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The lower bound of the shear viscosity to entropy density ratio is examined using an exact 
representation of the ratio through the density of states. It is shown that the lower bound in 
a generic physical system is not universal, its value is determined by the entropy density. Some 
examples of physical systems are discussed in the paper where one can expect violation of the 
conformal 1/47T value. 
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Experimental evidence from RHIC [1, 2] suggests [3, 4] that the QCD matter near the (would-be) critical temper- 
ature is very close to an ideal fluid, the shear viscosity to entropy density ratio is rj/s < 0.2 [5, 6]. Such small values 
J> ■ suggest very strong interactions, and it is therefore a big theoretical challenge to explain the actual i]/s in this system. 
Direct calculations with weak coupling resummed perturbation theory in various models [7-11] and resummations 
with Boltzmann-equations [12, 13] suffer generally from the problem that the unperturbed system, ie. the free gas 
has infinite viscosity, and so small perturbations necessarily remain in the large viscosity regime. This is reflected in 
the result rj/s ~ 1/ '(g 4 log<?) which is large if the coupling is small. At stronger coupling subleading corrections may 
alter this result: for example the inclusion of 2 — 3 Bremsstrahlung processes into the Boltzmann equation results in 
a significant difference near the critical region of QCD [14, 15]. 

Since perturbation theory is not well controlled in the strong coupling regime, one has to look for alternative 
methods. A potentially exact method is lattice Monte Carlo simulations. However the actual measurement of the 
transport coefficients is plagued by the necessity of analytic continuation from imaginary to real time based on a 
discrete set of data. The inversion of the integral equation for the analytic continuation then has very little sensitivity 
to the important infrared physics [16]. Therefore existing MC results [17] depend on not fully controlled assumptions. 

A completely different strategy is when one uses a dual description of the theory, where the strong coupling limit 
is mapped to the small coupling limit of the dual theory. In case of QCD-like theories this is possible only with a 
much larger symmetry group, the N = 4 SYM theories, where supersymmetry together with the holographic principle 
makes possible to work, at large JV C and large t'Hooft coupling, in a five dimensional gravity with AdS metric [18]. 
There the shear viscosity of the conformal field theory can be calculated from graviton absorption [19], resulting in 
C*~> ' the celebrated 1/47T. In a large class of AdS models this seems to be a lower bound [20, 21]. 

Although Af = 4 SYM is not QCD, one can argue independently that there is a lower bound in QCD itself [22]. 
— — ■ Then the value of this lower bound may come universally from the large N c large t'Hooft coupling supcrsymmetric 
0^ \ case [19]. 

The statement that 1/47T would be the lower bound in any system for the rj/s ratio is, however, not really proven, 
| only argued for, and it is supported by the fact that we never encountered systems where this bound was violated. 
But the argumentation has weak points, which suggest that it is not so general as it is commonly believed. First, 
in the gravity side one can find higher curvature AdS models, where this conclusion fails [23-25]. In the framework 
of nonrelativistic gases one can also construct counterexamples [26, 27]. Although these examples may describe 
mctastablc states [28] , but the question of the applicability range of the result is not answered reassuringly. 

More importantly, the argumentation is strongly based on the quasiparticle picture. The generic argumentation 
says [19] that rj/s ~ Et where E is the energy of an excitation, r is the lifetime, and uncertainty principle states 
that Et > h which yields rj/s > fi. But uncertainty principle, in fact, states that the uncertainty of the energy 
measurement is related to the lifetime, not the energy value itself. In quasiparticle picture the uncertainty cannot 
be larger than the energy value itself, but in a generic theory this may not be the case. So, in fact the restriction 
Et > Ti is not a law of nature, but rather the applicability range of the quasiparticle picture. If the small width 
quasiparticle approximation cannot be applied, then we have no real arguments in favor of the lower bound. Recently 
in [29] the authors have shown that the presence of a continuum in the particle spectrum has important effect on the 
shear viscosity to entropy density ratio. 

To clarify these issues one would need exact statements about the lower limit of rj/s in physical systems. This is 
what is attempted in this work, using first principles. For the derivation we have to make some assumptions about 
some current matrix elements, then one can show that if the transport is dominated by a finite number of quantum 
channels (ie. there are no large number of particle species) then at finite entropy density the rj/s ratio has a lower 
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bound. This lower bound, however, depends on the entropy itself; for small entropy density case it is 
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where Nq is the number of dominant quantum channels, L is the "interaction range" . This latter can be defined so 
that parts of the system separated more than L do not interact, or interact weakly; in other words in systems larger 
than L the entropy is an extensive quantity. For more details and for general s formulae see Section III and (46). 

The mathematical reason why rj min ~ s 2 depends qualitatively on the fact that, as it will be explained in the paper 
in detail, the transport coefficients come from integrals depending quadratically on the energy density of states g, 
while the entropy density, for low values of s, depends linearly on g. Since these integrals have positive kernels, we 
can treat them as some averaging. So, up to normalization factors 77 ~ (g 2 ) and s ~ (g). This implies that with this 
normalization i] > s 2 , and the minimum is therefore proportional to s 2 . 

In a finite degree of freedom quasiparticle system the entropy density is proportional to T 3 and the interactions 
are screened at most at a scale of the temperature, and so the right hand side is constant. The value of this constant 
cannot be determined from a general argumentation, but by the argumentation of [19, 22] it may be 1/47T. But as a 
mathematical statement we can claim that there is no lower bound for the shear viscosity to entropy density ratio. 

The above formula also tells us where can we expect to have small r//s ratio: if the number of dominant quantum 
channels is large (cf. [26, 27]), if the system interacts on very large scales (ie. large width systems, where the width is 
not connected to the temperature), or if we are at low temperatures and the entropy density vanishes faster than T 3 . 

In the paper we discuss some systems, where possibly violations of the l/4ir value are present. These are systems with 
strong wave function renormalization, large width finite temperature systems and low temperature non-quasiparticlc 
systems. All of these cases can appear in real systems: strong wave function renormalization shows up in strongly 
coupled quasiparticle systems; at finite temperature QCD we expect large width systems [30]. Finally in systems 
where there are zero mass excitations, the spectral function must not have a separated Dirac-delta peak even at zero 
temperature. 

The structure of the paper is the following. We first derive exact expressions for the Kubo formula and for the entropy 
(Section II) through the density of states. In Section III we study the mathematical structure of the 77/s ratio and 
determine its minimum. Next we study the quasiparticle-like systems and discuss the small-width quasiparticle case 
in Section IV. In Section V we examine different situations where the excitations are not small-width quasiparticlcs, 
and discuss the impact of the off-shell effects on transport. We finish with conclusions in Section VI. 

II. THE CURRENT-CURRENT CORRELATOR AND THE ENTROPY DENSITY 

In this section we derive exact representations for the Kubo formula and the entropy density using the density of 
states (or energy spectral functions). 

A. Kubo formula 

The shear viscosity, the linear response coefficient for the transversal momentum difference in a plasma, is defined 
through the Kubo formula in the framework of quantum field theory. We introduce 

C(x) = ([ir xx (x),n xx (0)]) , m = T tj - ^SyTf, *,j = 1,2,3 (2) 

where T is the energy-momentum tensor of the system, and the expectation value has to be taken in the equilibrium 
system: 

(Q) = ^Tre-P H Q, Z = Tre^ H . (3) 
The viscosity is obtained after Fourier transformation as 

Cbo,P = 0)_ 

This discussion can be generalized by using instead of 7Tjj a general conserved current Jj, and look at the linear 
response function 

Cj(x) = ([J i (x),J i (Q)]), (5) 
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(there is no summation on indices i). Because of rotational invariance of the equilibrium system, only the diagonal 
elements are nonzero here. The corresponding transport coefficient will be called r\j. 
We continue by inserting energy eigenstates into the expectation value above. 

Cj(x)= -| J2Un\e-^ H Ji(x)\m) (m\Ji(p)\n) - {m\e^ H J,(0)|n) (n\Ji(x)\m) 

n,m ^ 

= 1 E ( e ~ PEn - z~' iEm ) (n\Ji{x)\m) (m\M0)\n) . (6) 

The possible states in the system is characterized by their total energy k total momentum k and other quantum 
numbers JC 

\n) = \k ,K!C). (7) 

Similarly we denote \m) — |<7o>q> 2)- We will also use the four-momentum notation k = (fc ,k) and q = (qo 7 q). 

In a realistic quantum field theory the allowed energy eigenvalues can form a continuum, and at finite temperature 
this is always the case. Therefore it is more adequate to write the sum over the energy in the above expression as 
an integral over a spectral density. In finite volume spatial momenta form a discrete set, but at larger volumes a 
continuum representation is applicable for them, while other quantum numbers remain discrete: 

/dkc, t d^k 

\k,K){k,K\ = vY J /7^I^( fc ) \k,lC)(k,lC\. (8) 

Since the state \k, JC) is dimensionlcss in our normalization, therefore g has the dimension of inverse energy. With the 
energy density of states we can write 

^W = ^E/(04^( fc )E/(04^W {e- pko -e-^)(k,K\J i (x)\q,Q)(q,Q\J i (0)\k,lC). (9) 

Using the generator of the space-time translation we write 

(k^JiWlq, Q) = e< k -^ x (k,K\Ji(0)\q, Q) , (10) 
then the Fourier transformation yields 

^(p) = ^E/70^( fc )E/7^y4^W(2^ (4) (fc-? + P) (e-^-e-^)\(k,K\Ji(0)\q,Q)\ 2 . (11) 

Then the transport coefficient reads 

^ = /^E l^^(k)0Q(k)e- f3ko \(k,JC\J t \k,Q)\ 2 . (12) 

Since Ji cannot change the quantum channel without making any change in both energy and momentum, so nonzero 
result comes only from JC = Q: 

^ = ^Y t j^^)e- pko \{k,K\Ji%K)\ i . ( 13 ) 
The (temperature-independent) matrix element of Ji can be calculated in the free theory with the result 

(k\Ji\k) = q-^ (14) 



where q is the charge carried by the current. In case of the energy momentum tensor q = k.y. 

Vk 



(k\T xz \k} = ^. (15) 
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In general, we expect an expression of similar structure, since the external momentum carries the Lorentz index 
of the current, and by dimensional reasons it must be divided by an energy- like quantity. There can still be a 
Lorentz-invariant factor: 

(k,IC\Ji\k,IC}=q 1 ^J K (k 2 ), or (k, IC\T xz \k, AC) = ^ T K (k 2 ). (16) 

The reduced current matrix clement J7x; depends on the system under consideration, and we cannot tell much more 
on a generic ground. We will assume, however, that Jk. does not grow at large fco exponentially. 

Since the only directional dependence is in the current, we can average angular dependence using (16). Since 

cos 2 = 1/3 and sin 2 6 cos 2 9 cos 2 <f> = 1/15, wc can write 

"-s «*<*>«<*»'■ (17) 

and for the shear viscosity 

'-izZl&^W <*<*>"=<*»• < 18 » 

The volume dependence canceled, and so we interpret these formulae as infinite volume limit results. 



B. The entropy 

To calculate the entropy density we start by computing the free energy as 

Z = e-0* = Tre"^ = £ = ]T / ^(fc)^ = V QK{k)e~^. (19) 

In a strongly coupled system the volume dependence of the free energy can be arbitrary. When we increase the 
linear size of the system, after a certain scale L one can observe more and more accurately the the linear volume 
dependence, ie. the free energy becomes an extensive quantity. That means that adjacent volume elements larger 
than L interact (mostly) with surface interaction, while smaller volume elements interact through the total volume. 
Therefore this "interaction range" serves also as an infrared cutoff for the relevant interactions. In quasiparticle 
systems the interaction range corresponds to the linear size of the cross section, the inverse of free mean path. In 
general strongly interacting systems have large L, weakly interacting ones have small L. 

This is the scale, beyond which a coarse grained description is sensible [31, 32], since only there can one speak 
about densities of different physical quantities. The free energy in a coarse grained description is dominated by the 
integral of the free energy density, while the corrections coming from the surface interaction are subleading. Therefore 
in a microscopic description we should choose the volume size at least V = L 3 in order to catch all the relevant local 
physics effects. Therefore we define our free energy density as 

/--^'E/lSF «««-*■)• (20) 

We can separate the ground state (vacuum) contribution, where always [40] g vac {ko) = 2irS(ko): 

/.^(l + L-E/^W.-*.). (21, 
From this generic form the entropy density reads 

s = TE ln ( 1 + i3 E/(0 acOfe) e- 0ko ) +|E/(0i^o acOfe) (22) 

\ JC / K. 
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Then the transport coefficient to entropy density ratio reads 



In case of shear viscosity we have to substitute J 2 — > k 2 T 2 /5. 



(23) 



III. MINIMIZATION OF SHEAR VISCOSITY TO ENTROPY DENSITY RATIO 

Equation (23) gives the functional dependence of rjj/s on £>, and so we can study mathematically the minimum 
of this ratio by varying g. The above formula, however, contains yet the unknown J factor. To say something 
constructive we have to make assumptions on this function: we will assume that it can get at most power law ko 
contributions in a strongly interacting system. In this case we can perform an analysis to find the lower bound. 

There is a simple way to decrease the value of r/j / s: wc start with an arbitrary g with finite r/j / s and make a 
rescaling g — > Xg with A — > 0. Then the integrals of g become small, we can make a linearization in s with respect to 
g. Then the numerator is quadratic, the denominator is linear, and so rjj/s — ► Xr/j/s, in the A — > limit it vanishes. 

In a physical system this simple procedure cannot be performed, because we have a sum rule in each quantum 
channel 



/ 



^ftc(*0,k) = lfc(k), VAC.k. (24) 



But this constraint is not really restricting. To understand it we should remark that all g integrals for rjj and s are 
weighted by e~^ ka . Therefore rjj/s is sensitive to the ko < T regime, say, the infrared (IR) regime. The sum rule, on 
the other hand, depends also on ko > T regime, the ultraviolet (UV) regime. Therefore we may rescale the IR part 
of the energy density of states and still maintain the sum rules. 

To sec this we assume that we can linearize the entropy density in g and try to find the minimal r\j. The generic 
case we will discuss later. The assumption of linearizability also means that Z w 1. To simplify the notation we 
introduce 



E 



d 3 k 



(25) 



(2tt) 3 

For the linearized case, in this symbolic notation the transport coefficient and the linearized entropy density read 

VJ = J J^a Q (k )g 2 Q (k )e-^, s = J J ^ 7^0)^0)^°, (26) 



IQJ Z7r JQ 

where 

/3k 2 
to 



a Q (ko) = ' w J£(k), 7 (fc ) = 1 + 0ko. (27) 



We shall minimize rjj with the constraint (24), and so we apply Lagrange multipliers to each quantum channel. To 
have a control to the thermodynamics we also fix the entropy density s. Then we have to minimize the expression 



mm 



'JqJ^ [^Q(ko)g 2 Q (k )e-^-2X 1 (k )g(ko)e-^ -2X Q g Q (k Q )]. (28) 
Its solution provides the spectral function 

A 7 (fc ) + Aqe^Q 

gQ{k0) = MM • (29) 



The normalization condition can be satisfied if 



_U R -XX 1 ,_ s -Iq UqXi/Xq 

Xq -^^> x - j Q (x 2 -x!/x y (30) 
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where 

e" ko v _ [dk 7 (fco) v _ [dk 7 2 (fc 



_ [dkp _ fdkp 7 (fco) y- _ [*rW ,„ n 

A ° 7 2^ a Q (k o y Al y 2^ ag(fco)' 2 y 2tt a Q (fc ) ' [ j 

Then the minimal ?7j reads 



x 2 - f] Jq Xo 

The concrete expressions for Xf are potentially divergent. To define still the result we have to introduce an energy 
cutoff A into the system. Then the qualitative behavior of the quantities, assuming J 2 {ko) <~ k$: 

X ~ Ag-°e^ A °, Xi-A 4 " a , X 2 - convergent. (33) 

Therefore removing the energy cutoff results in vanishing X" /Xo factors. This is exactly the mathematical appearance 
of the qualitative analysis in the beginning of the subsection: the sum rule constraints of (24) can be easily satisfied 
by tuning only the UV part of g. 

When we remove the UV cutoff what remains is 

s 2 

VJmin = ~r^r- (34) 

h X2 

Evaluating the integral in the denominator, for concreteness for shear viscosity where we assume J ~ k 2 , we obtain 

(35) 



NqLT a ' 



where Nq is the effective number of quantum channels. This formula is valid for small s. 

We remark here that the same line of thought leads to the conclusion that also rj/s has a lower limit which is 
proportional to e itself. Therefore s has no singled out role in this analysis. 

We can perform the same analysis for large s, too. This time we omit the sum rule constraints, but we cannot omit 
Z. Our strategy will be to fix Zs and minimize Zrjj. Therefore we have to minimize the expression 



/ J d ^^ Q {ko)el{k )e-^-2x(^+ j J^pk ogQ (ko)e-^ 



(36) 



where 



2 = l + L 3 |J^ Q (Me-*. (37) 

The minimum condition yields 

l + hiZ+j3ko 

Qdfa) = A — rpr . (38) 

«q(ko) 

If \nZ <C 1 we get back (29), this corresponds to the small s case. Here we follow the opposite limit and assume that 
Z dominates the above expression in the relevant k n <T domain: 



9 Q (ko) « A^^. (39) 
aQ(ko) 



With this density of states we can calculate 



Z = 1 + \Y \nZ, Zs = (Z\nZ + \Yj_\nZ), (40) 
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where 



Yq = L / / ^ 



dfc e-^ feo 



2?r a Q {k )' 



Yl =L 3 ^ 



Q 



2tt ag(fco) 



(41) 



For the case of the shear viscosity, assuming Jq ~ k 2 we obtain Y\ <~ Yo ~ Nq(LT) 4 . For J7q = k 2 there is a factor 
of 3 between them, so we may assume that their values are of similar order of magnitude in general. Then in the 
expression of Zs the first term is proportional to (InZ) 2 , the second only to InZ. Therefore the second term can be 
omitted, and we find 



Zs. 



Z\nZ 



s = 



InZ 



(42) 



Therefore in the large Zs regime s and Z are not independent. So if we fix a large value for Zs it means fixing large 
value for separately s and Z. Large Z also means that the vacuum contribution can be omitted, and then 



A 



Z 



Yo InZ 



Finally for the transport coefficient we obtain 



therefore the r//s ratio reads 



Z 



D L J s 



1 



in Y L 3 s L 3 s Nq(LT) 4 ' 
So we can summarize the small and large s cases in a common formula 



.F(L 3 s) 
N Q (LT) 



4> 



(43) 



(44) 



(45) 



(46) 



where for small values F(x) ~ x, for large values F{x) <~ e x /x. 

The results has several interesting consequences. First of all we see that if we distribute the entropy uniformly into 
Nq quantum channels, we can reduce the ij/s ratio at any temperature to zero. This corresponds to the construction 
of [26] with large number of particle species. In that model L 3 s ~ In Nq because of the mixing entropy, which for 
large Nq leads to 1/lniVg dependence for the shear viscosity to entropy ratio. In this construction the entropy goes 
to infinity, while the shear viscosity stays constant. 

If we fix the number of dominant quantum channels then the r]/s ratio has a minimum for each fixed entropy, but 
this minimum depends on the value of the entropy itself. The theoretical lower bound is zero, which, however, can 
be reached only at zero entropy, ie. at zero temperature. 

If we have non-interacting quasiparticles, then the effective interaction range is zero, L — > and the r\j s ratio is 
infinity. In case of weak coupling the interaction range is related to the cross section, or the inverse lifetime of the 
particle. A more detailed description for the quasiparticle systems is done in the next section. 



IV. QUASIPARTICLE SYSTEMS 

Apart from the theoretical bound determined above, it is still a question, how the r\j ' s ratio behaves in a real system. 
In this Section we actualize the general formulae for quasiparticle systems, and discuss the small width case. 

For quasiparticles we use that the energy contribution from non-interacting subsystems simply added in the free 
energy, so in weakly interacting bosonic/fcrmionic system we can write 

f=T Jw)~ i QQp[k) (T) ln (1 ± e ^ fe0) ■ (47) 

The entropy density can be calculated as s = —df/dT. 

QQ(k) = ^J-0^X(Pk o ) QQ (k). (48) 



E 



(27T)* 



-3^-±In(l±e-^) 

2 /3fe ± 1 V > 
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For concreteness we will calculate the shear viscosity coefficient from now on. The viscosity to entropy density ratio 
reads 

, is ?/<&^ ■-*<*<*>"<*»* 

1 IZ ■ (*5> 



Q 



Formally it is similar to the perturbation theory motivated expressions of [30]. But here the logic is different, we do 
not use propagating states, but instead energy eigenstates as intermediate states. 

To treat the problem analytically we will assume that J7jt(fc) does not depend too strongly on the momentum on 
the relevant momentum regime, where, for example g(k) exhibits a peak. We also will assume that there is a single 
dominant quantum channel where the largest contribution for entropy as well as the shear viscosity comes. Therefore 
we analyze a reduced shear viscosity from the dominant quantum channel as 

V - 15 J (2tt)* kl C Q ^ k >- (50) 

The rj/s ratio then reads 

P f d*k (k 2 ) 2 0kn 2 



S 



We hope that r\ characterizes well the true r]/s, so we will use this quantity to assess the importance of the off-shell 
effects. 



A. Small width approximation 



If the integral is dominated by a Dirac-delta-like peak then near the peak region we may approximate 

2r 



QK(k) 



(52) 



In relativistically invariant systems e\ — k 2 + m 2 , and approximately the same is true at finite temperatures with 
thermal masses. 

The product of two functions can be approximated by 



£»/cO) ~ £/c(fco = £k)g(k) w - 27T(5(fco - e k ). 



(53) 



By replacing back this expression into (50) we find 



f 



VQP 



15ttTT 



dk^e~^ T , 

4 



(54) 



where the QP subscript refers to the small-width quasiparticle case. Depending on the dispersion relation we can 
obtain results like 



VQP 



( 8 T 4 

1 m3/ 2 T 5 / 2 



-m/T 



2ir T 

4 m 7 / 2 T 1 / 2 

I i5V2^ r 



if Sk = k 



if s k = m + - 



2m 



(55) 
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The first case corresponds the high temperature case, the second the low temperature relativistic case, the third a 
low temperature nonrelativistic gas case (Chapman-Enskog formula). 

The free energy in the quasiparticle case (52) reads after partial integration: 



/qp = -^2 JdEk(Efn(E), 

Eq 

where n±(E) = (e@ E ± l) -1 Bose-Einstein or Fermi-Dirac particle number density. For the entropy we obtain 



(56) 



OO 



(57) 



In the different dispersion relation case we get: 



sqp 



E 



4T 3 TT 2 



A± 90 

yl/2^5/2 

(2tt) 3 / 2 ' 
{mTf 2 
l2(27r) 3 / 2 ' 



-(3m 



if £fc = k 



if £ k = m + £j , 



(58) 



if Sk = k 



2m • 



where A± = (1,8/7). Then the ratio reads 



sqp 



36 T 



= < 2tt 



a±tt 4 r ' 

16?r m 2 
~15 FT' 



if £fc = k 



(59) 



In the first, m <C T case the width must determined by the temperature which is the only scale in the system: 
r ~ T by dimensional reasons. Therefore in the first case we obtain temperature independent ratio. 

In the second, low temperature massive case we have to take into account that the quasiparticles are stable at 
zero temperature, and the first scattering states which are responsible for the decay have to be excited thermally. 
Therefore we expect that T ~ e~ M / T , where M is the mass of the excited state. In this case the ratio grows like 
T 2 e M / T at small temperature. 

In the third case the width can be proportional to T, but also in this case the ratio will grow like 1/T 2 . 

So we can see that in the quasiparticle case the rj/s ratio is at most constant, but in all massive cases it is bounded 
from below. The lower bound of this ratio is probably provided by the AdS/CFT results. 

The inverse proportionality to T also suggests that in environments where the correlation lengths, in particular the 
quasiparticle lifetimes grow, the i]/s ratio will decrease. This is the case for second order phase transitions. 



V. RELEVANCE OF THE OFF-SHELL EFFECTS IN DIFFERENT PHYSICAL SITUATIONS 

In the field theoretical systems besides the (quasi)particles one should count with a multiparticle contribution which 
appears as a continuum in the energy density of states. The simple picture suggested by the small width quasiparticle 
approximation of the previous subsection may break down more or less in case of strong continuum. 

In this section we consider three physical situations where the presence of off-shell effects can significantly modify 
the small width quasiparticle picture. The first is the effect of the wave function renormalization, the second is the 
case of large thermal width, the third is the case of a theory with zero mass excitations. 



A. Wave function renormalization 



As a first example we take a system where besides the quasiparticle peak there is a separate continuum. The 
normalization requirement (24) tells us that in the full case the peak should be smaller than in the quasiparticle 
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approximation: the factor we have to apply is the wave function renormalization. At small enough temperatures we 
will see only the quasiparticle peak, then the rj/s ratio is simply rescaled by the wave function renormalization factor. 
To demonstrate the effect we take a toy model, where the full spectral density is taken as 

g(k) = Zg QP {k) + 2tt 1 - - 0(E 2 >k > E{), (60) 
tLii — tb\ 

ie. we take a step-function-like spectral density in the continuum regime, which means that in this toy model we do 
not care about the details of the continuum, only an effective height is taken into account. Here we take the relativistic 
form E\ 2 = k 2 + M 2 2 , while qqp has the form of (52). The proportionality factors are taken into account to satisfy 
the sum rule 

/f%(M = /|%QP(M = 1- (61) 

We will assume the hierarchy m <§C M 12 and T < M\. In this limit in the square g 2 we can neglect the cross terms, 
and we can use nonrelativistic approximation in the £1,2 energies. What we obtain in the leading order 

where AM = M2 — M\. With similar assumptions the entropy density reads: 

1 (1 - Z)T 3 / 2 Mf /2 M /T 
3 = ZS * P + (2^- AM e . (63) 

If the temperature is lower than the threshold, then the radiative corrections for both fj and s are suppressed. Then 
approximately we obtain 

^=Z^. (64) 
s s QP 

That means that even though the continuum contribution cannot be seen, its effect is measurable in the rj/s ratio. 

The value of Z is determined by the relative contribution of the continuum in the sum rule as compared to the 
quasiparticle peak. In weak coupling cases this ratio is suppressed by (powers of) the coupling constant, so Z w 1 and 
so the effect on rj/s is small. However, in strongly coupled theories the continuum becomes more and more important, 
and, correspondingly, the wave function decreases. As Z — > the continuum contribution slowly takes over and it 
will dominate the rj/s ratio. In this simple case the theoretical Z — > limit the continuum yields fj/s ~ T/ AM. This 
demonstrates a mathematical example where in a constructed system the rj/s ratio can go to zero. Still, we have 
to remark that there does not exist a physical system where this type of continuum would be manifested. But it 
emphasizes the importance of the continuum part of the spectral function in the transport. Two more realistic cases 
will be studied in the next subsections. 



B. Large width case 



In lots of realistic examples there is a quasiparticle peak, but its width is large, comparable to the mass scales of 
the system, cf. for example [30]. At finite temperature strongly coupled theories, moreover, the quasiparticle peak 
can merge with the continuum, forming a broad, slowly varying spectral function, as it can be seen for example in 
the 2PI simulations [33-35]. In this case the "width" of the peak has nothing to do with the quasiparticlcs. In fact 
there are no real particles in the plasma, the excitations decay before a particle is formed. Correspondingly, as the 
width grows, the quasiparticle picture becomes less and less good approximation. 

To treat this effect numerically in the case of the transport coefficients, we take a spectral function which could 
describe a spectral function of a quantum channel at finite temperature, and which is of a Breit-Wigner-type form 

eBw[k) = * + if (65) 



The Af prefactor is necessary for normalization, is Af — (tt + 2 arctan(_E 2 /7 2 ))/(27r). In the above formula the 
dispersion relation and the momentum dependence of the 7 parameter is arbitrary, but we assume spatial rotational 
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invariance. In the weak coupling limit the width corresponds to T — 7 2 / (2Ek). In the high temperature field theories 
we expect 7 ~ g 2 T. In strongly coupled case 7 3> T is also possible. 

In the fj/s formula we rescale by the temperature, so we choose T = 1. Then we have to calculate 



1 

s 



dkk b J dk Q e~ k ° (QBw/koY 
i' 



(66) 



15 Jdkk 2 J dk (fc (e ° - I) -1 - hi(l - e- ko )) g BW 


where we have chosen, for the sake of concreteness, the minus sign in X . 

For a particular choice of 7*. = 7 = const and = ^/k 2 + m 2 we find the plot of the ratio in Fig. 1. This figure 
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FIG. 1: The shear viscosity to entropy ratio as a function of the imaginary self energy. 



shows how the small width quasiparticle picture which yields in this case ~ T 2 /^ 2 breaks down for larger width 
case, where we find ~ j 2 /T 2 behavior. Looking at these figures as function of temperature at fixed 7 we see that 
at high temperatures the r\j s ratio increases with increasing temperature, while at low temperature it decreases with 
increasing temperature. While the former is characteristic for gases, the latter is the behavior of the fluids. So in this 
approximation we can give an account for the 77/s ratio in the fluid-gas crossover. 

How can we understand mathematically this behavior? The small width regime corresponds to the quasiparticle 
approximation above, so we have to understand the large 7 regime, so examine the 7 » T case. Because of the 
exponential suppression for ko, the value of fco cannot grow larger than T. The peak region, however, is in the vicinity 
of |fcg — E 2 \ ~ 7 2 . Therefore fco cannot play a role in saturating the integrand (unless we are in the quasiparticle 
regime 7 -C m), and so we can neglect kg here. Then the dominant contribution should come from Ek ~ 7, since 
there is no exponential suppression for k. As a result fco <C k and the fco and k integrations decouple. For )»ra 
case fit ~ 7 > m means that we can use ultrarelativistic dispersion relation. Then for the numerator we have 



dko e 



dk k b 



1 



I67 4 



A/-2 Qgi + 7 4)2 



(67) 



and for the denominator we obtain 



dkok (k (e ko -I)- 1 -\n(l-e- ko )) dkk 



fc 

Therefore we expect that the ratio grows like 7 2 /T 2 . In the opposite case, when 7 <C m, but the decoupling of the 
integrals is still true, then we can neglect the 7 factor in the denominator and we obtain 



dkk" 



1 



I67 4 



7 



J\f 2 {{k 2 



(69) 
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and 

oo oo 

J dk Q fc (fc (e fc » - I)' 1 - ln(l - e- k °)) J dk fc 2 l (fc2 ^_ „ (70 ) 



Therefore the ratio behaves as j 2 /m, so the trend of large 7 regime continues. 

To have a feeling of when does the small-width description start to dominate, we remark that the quasiparticlc 
contribution is at one hand suppressed by e~ Ek , but, on the other hand, the denominator can be small, too. Decoupling 
of the integrals will not be true if the denominator can compensate the exponential suppression. In the small mass 
case, when also 7 » m is true, the quasiparticle regime gives e^ k /'j 8 with k <~ 7, as opposed to the I/47 8 factorized 
case. So on-shell effects will become strong if 7 < 1.4. In the large mass case the on-shell contribution is e~ Ek /j 8 
with Efr <~ m, the off-shell effects yield 1/m 8 . Then the borderline is at 7 <~ me"" 1 ' 8 . In the hypothetical m — > 00 
case, therefore, there is no quasiparticle regime, the linear 7 2 trend will continue until 7 = 0, where the ratio is zero. 



C. Low temperature systems with zero mass excitation 



If we wonder whether we can reach in a real system zero value for 77/s, with finite number of quantum channels, we 
must reach, according to (46), zero entropy. In a real system can happen only at zero temperature. In the small- width 
case, as we have seen, the 77/s ratio diverges in the massive case, and, formally, stay constant in the zero mass case. 

The zero mass systems, however, need a more detailed study, since there are never classical Dirac-dclta-like spectral 
functions. To see this, we recall that if we take a quantum field theoretical system with a stable particle with mass m 
which interacts with another particle with mass M then the spectral function of this system will contain a Dirac-delta 
peak at kg = m and a continuum starting at the threshold m + M . The imaginary part of the self energy, because of 
kinematical reasons, is proportional to yj A(p 2 , m 2 , M 2 ), where X(x,y, z) = (x — y — z) 2 — Ay 2 z 2 . For a finite M this 
leads to a square-root behavior for the complete spectral function, g ~ \Jp 2 — (m + M) 2 near the threshold. 

If the massive particle interacts with a zero mass par ticle, all this means that there is no gap between the quasipar- 
ticle peak and the continuum. Moreover, in this case = p 2 — m 2 , which is the same as the tree level part 
of the real part. This has the consequence that near the threshold the spectral function is divergent g ~ (p 2 — m 2 ) -1 . 
There are logarithmic corrections to this behavior; for example in QED at one loop level in Lorenz gauge [36] 



g(p 2 k, m 2 ) 



2Q(p 2 -m 2 )- 

7T 



(j/> — m) 





2 2 
p — m 






) 




M 2 





(71) 



The leading logarithmic corrections can be resummed by renormalization group or by Bloch-Nordsieck construction 
[37], the result is the modification of the threshold behavior 

. 9 9 , 20(p 2 — m 2 )(i> + m)m~ 2 P sin a a e 2 

The main features, however, still survived: there is no gap between the quasiparticle peak and the continuum, and 
there is a rising spectral function as we approach the mass shell. 

This behavior can be studied numerically, too. In case of QCD one can use Monte Carlo data to reconstruct the 
spectral function [38], which exhibits qualitatively the same behavior as in the QED case. 

According to the above analysis in all systems containing zero mass particles, the energy density of states cannot 
contain a zero width quasiparticle. The would-be quasiparticle continuously radiates long-living soft zero mass exci- 
tations, and continuously interacts with them. The spectral function is not Lorentzian, it has no "width" in the usual 
sense, and it does not describe a decay, too. At larger temperature we will still observe a quasiparticlc, although 
the mechanism how it is formed is far from trivial [39]. At low temperatures, however, we are in a non-quasiparticlc 
system, where the effect of the threshold behavior cannot be neglected. 

At small temperatures we are in an almost Lorentz invariant system. Because of the volume normalization we can 
write g{k) = k g(k 2 ), where g depends only on the Lorentz invariant k 2 form. Then we write (cf. (51)) 



15 j (^^(/SfcoK^ ^) 



(73) 
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The exponential factor forces the system to have the possible lowest energy values, ie. near the threshold. 

First let us assume that the lowest lying threshold is at k 2 = M 2 with finite mass. This is conceivable if the 
zero mass particle is not asymptotic state, like in low energy QCD. Near the colorless bound states, which have 
color multipole moments, the gluons can still exist, and the effective gluon cloud can result in a continuous threshold 
behavior. In this case we can use (3k > (3M » 1, and so X(f3k ) ~ (3k . We change into 4D polar-coordinates 
(ko = fc cosh 77, k z — fcsinhrjcosfl, k y = k sinh 77 sin sin <j>, k x = k sinh 77 sin 8 cos <fi) , taking into account that ko > 0, 
and the 3D rotational invariance: 



d4fc rv/, ^ 1 



(2tt)' 



dkk 3 / di] sinh 77. 



(74) 



We use the integral formula 



/ 



d? ? sinh r V-* cosh7 ' 



V^(z/2) n / 2 



(75) 



(K n is the modified Bessel function of the second kind), and the asymptotic form K n (x) — yjir/{2x)e x to arrive at: 



T 2 J dk k 7 / 2 e~ l3k g 2 (k) 

V = __o 

„ 00 

/ dkkV 2 e-P k g{k) 



Near the threshold we power expand the spectral function like 

g(k ?a M) — CQ{k - M) {k 2 - M 2 ) w . 



(76) 



And then find 



fj _ r(i + 2w) 
7 _ r(i + w) 



C(2M) w T 2+y 



(77) 



(78) 







This result cleanly shows that if there exist a system described above, it must have a vanishing rj/s ratio in the T 
limit. 

If the zero mass particle is an asymptotic state, then the lowest lying energy eigenvalues belong to the quantum 
channel of the massless particle. Since it interacts with itself, the spectral function has no gap here, too. A special 
class is the conform theories, another the weakly interacting gauge bosons, like the photon gas, where the photon- 
photon scattering is mediated by virtual electron loop. In the massless case the above analysis goes through without 
modification and it yields 



- - CT 2{1+W) . 
s 



(79) 



This result could be obtained by purely dimensional analysis, since the dimension of C is —2(1 + w). The numerator 
of (73) contains g 2 , the denominator only g, so there remains a factor C in the ratio. Since the ratio is dimcnsionlcss, 
something has to compensate this factor. In the massless case, only the temperature can do that: thus we obtain 
~ T 2 ( 1+w ) dependence. 

In conformal field theories without anomalous dimensions C is dimensionless and w = — 1. That predicts a 
temperature-independent r\j s ratio. In other cases, ie. either in a non-conformal field theory or in a conformal 
field theory with anomalous dimension w can differ from —1, then we again observe a vanishing 77/s ratio in the zero 
temperature limit. 



VI. CONCLUSION 



In this paper we examined the shear viscosity to entropy density ratio using exact representation through the 
density of states or energy spectral functions. We examined what can be said purely mathematically about the ratio, 
assuming some physical conditions (sum rules) for the spectral functions, and keeping the entropy density constant. 
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We concluded that the 77/s ratio has no lower bound in a most generic class of physical systems. To understand this 
statement qualitatively we recall that r\ ~ g 2 and s ~ g if the temperature is low enough. Therefore if the density of 
states exhibits large peaks (quasiparticle systems) then g 2 3> g and so the rj/s ratio is large; if g is small everywhere 
then g 2 <C g and the ratio is small. However, if we fix the entropy density then there is a minimum, where r\ <~ s 2 . 
This means that rj/s can go to zero (in a fixed system) only at zero entropy density, ie. at zero temperature. 

Although mathematically the situation is clear, it is hard to show physical systems where a small r\j s ratio can occur. 
First of all if it consists of quasiparticle at all, they must not have small width, since the small width approximation 
excludes small rj/s ratio. Candidates for such a system are quasiparticle systems with strong continuum (small wave 
function renormalization constant), high temperature strongly interacting systems or low temperature systems with 
zero mass excitations. In the first two cases it is possible to go below 1/4-7T, but we can never reach zero value. In the 
third case if the threshold behavior is (k 2 — M 2 ) w then the r\j s ratio can reach zero as 7 l2 ( 1 + u ') ( m M = case) or 
T 2+w (massive case). So if w > — 1 then we will find a vanishing ratio in the T — > limit. 

As far as the heavy ion experiments are concerned, the QCD at finite temperature and density represents a finite 
entropy density system. Therefore there is a lower bound for the rj/s ratio, but its value is not necessarily 1/47T. 
However, the relevant scales are set either by Aqcd or the temperature (since the couplings are of order one), and 
these two scales are again similar. Therefore we are close to a one-scale system, where the predictions of the conformal 
theory may apply. 
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